Journal of Molecular Evolution

Gene duplication is an important evolutionary mechanism contributing to genomic innovation and adaptation. While the long term fate of duplicated genes has been extensively studied, immediate effects such as those resulting from increase in gene dosage remain understudied. Of particular interest is the question of how expression level relates to the optimum after a duplication, as it could fall short, exceed or directly reach the optimal expression. In this study, we investigated the relationship between expression and fitness for 11 genes previously reported to increase fitness when duplicated in Saccharomyces cerevisiae to compare the expression level achieved by duplication to the one maximizing fitness. While most of the studied genes show little fitness variation across the expression levels tested, we found that ESP1, which encodes a separase essential to anaphase initiation, exhibits a clear fitness optimum at an expression level above the native single-copy state. However, the expression level associated with gene duplication exceeds the optimum. We find that ESP1 overexpression is associated with a reduced cell size and a lower proportion of the cell cycle spent in G2 phase, consistent with its role in cell cycle progression. We propose that the duplication of ESP1 accelerates the cell cycle, resulting in smaller, faster-dividing cells. When measured in different environments, the fitness functions also exhibit similar shapes, with slight variation in the optimal expression. Our findings highlight how gene duplication can confer an immediate selective advantage through dosage increase while being non-optimal, reversing the directionality of natural selection on expression level after duplication.

Single-cell sequencing provides a rich source of data that can be used to model and estimate the evolution of somatic cells within individuals. Many common models for modeling DNA evolution in phylogenetics are based on continuous-time Markov processes that operate on aligned sites. For haploid DNA, such as found in mitochondria, there are many well established models for these processes but since nuclear DNA is diploid, containing copies of both the maternal and paternal genome, it is natural to develop continuous-time Markov processes adapted to diploid DNA. One possible model is that both the maternal and paternal evolve independently according to haploid models so that, in principle, DNA comes with a phase. i.e. a label that tells if the site was derived from the maternal or paternal sequence. However, this information may not be available in practice and only an unphased diploid genotype is available. That is, an unphased model groups both of the possible maternal-paternal labels into a single state. This grouping process, called ``lumping'', is used by other methods as well and results in a coarser process. Care must be taken so that the coarser process can still be modeled as a continuous-time Markov chain; the lumped process of the original phased Markov chain may or may not be Markovian, depending on restrictions imposed on the generator matrix. Here, we explore the mathematical formulations of phased and unphased Markov processes, as well as a collection of other diploid DNA models from the literature.

Cells and tissues have the ability to respond to different stresses and maintain their homeostasis. The b-Zip transcription factor Xrp1 was found to play a pivotal role in this context in Drosophila. Xrp1 is required for the response to irradiation, in genetic conditions of ribosomal protein (Rp) deficiency leading to cell competition, and more generally in mediating the Integrated Stress Response (ISR). However, how these diverse stresses converge to and activate Xrp1 is still unclear. Our results demonstrate the existence of three separate mechanisms of Xrp1 induction. In response to Rp stress, a short isoform of Xrp1 is produced and requires the function of the ribosomal protein RpS12. In response to irradiation, a specific mRNA isoform is transcriptionally induced, which encodes the Xrp1-short protein isoform, in a p53-dependent manner. In response to the ISR, the long isoform of Xrp1 is translated by the derepression of a uORF present on the Xrp1 5'UTR, a mechanism that requires eIF2 phosphorylation, RpS12, and a transcriptional induction by ATF4. Therefore, our work provides a comprehensive understanding of the role of Xrp1 in different situations of cellular stress, illustrating the complex level of integration of the mechanisms controlling this central gene.

Two rounds of whole-genome duplication (2R-WGD) in the vertebrate lineage are thought to have driven gene diversification. Mammals have four genes encoding five nebulin superfamily proteins (LASP1, LASP2/NEBL, NEB, and NRAP) hypothesised to have arisen from invertebrate LASP via WGD. Here, we examined the evolutionary fate of the nebulin superfamily after an additional WGD, teleost-specific genome duplication (TGD). Genome and transcriptome analysis of zebrafish (Danio rerio) and medaka (Oryzias latipes), of which the complete genomes were recently determined, showed only lasp family genes retained duplicate copies. Zebrafish retained two LASP2/NEBL-derived ohnologs; one (LASP2L) had been misidentified as a nebulette gene in the Zebrafish Information Network (ZFIN) despite lacking nebulette-specific exons. Another (LASP2/NEBL) is newly identified in this study because it lies in a newly determined region of the latest genome assembly. Contrastingly, medaka retained two LASP1 ohnologs and one LASP2/NEBL. The duplicated lasp family proteins showed overlapping and specific tissue expression patterns. We confirmed all three zebrafish lasp proteins bound actin and co-localized with actin-rich structures when ectopically expressed in cultured fibroblasts. This supports that TGD-derived lasp duplicates underwent lineage-specific retention and expression divergence in teleosts whilst conserving core actin-binding properties.

Forms I&prime, I&Prime, and I RuBisCOs were recently identified in bacterial metagenome-assembled genomes (MAGs) and are significant for evolutionarily bridging form I, which is responsible for the major part of carbon fixation on Earth, and other forms. Pioneering studies focused mainly on the biochemical and structural characteristics of "form I-related" RuBisCOs and deduced the evolution of this protein family over phylogenetic relationships among form I and form I-related sequences. We here surveyed form I-related sequences in public databases and conducted phylogenetic analyses incorporating the updated diversity of form I-related RuBisCOs. This study yielded two major outcomes regarding RuBisCO evolution. First, the analyses presented here suggest uncertainties overlooked in the previously published studies. It appeared to be difficult to conclude (i) the phylogenetic relationship among forms I, I&prime, and I&Prime, (ii) the monophyly of form I&prime, or (iii) the monophyly of form I. Second, we identified a gene encoding form I RuBisCO in an archaeal MAG, and this protein is likely involved in the pentose bisphosphate pathway. The results presented here suggest that the evolutionary route leading to form I in RuBisCO evolution remains unsettled.

Prospero-related homeobox (Prox) genes encode a family of transcription factors that play essential roles in the development of several organs and systems, including the central nervous system, lymphatic endothelium, musculature, and liver. Despite their developmental importance, the evolutionary history of the vertebrate Prox gene family remains poorly understood. In this study we combined phylogenetic and synteny analysis to characterise the evolution of the Prox family in vertebrates. Our results reveal that two to three Prox subfamilies were already present in the last common ancestor of jawed vertebrates. We clarify the identity and evolutionary relationships of well-studied members of this family and identify multiple independent losses of Prox2 and Prox3 genes in specific vertebrate lineages. Furthermore, we uncover evidence for the existence of a fourth Prox gene in the ancestral vertebrate genome, which was subsequently lost. Overall, this study provides the first comprehensive analysis of the evolutionary history of the vertebrate Prox gene family and establishes a foundations for future studies on the functional roles of these genes.

The teleost whole-genome duplication (TGD) contributed to functional diversification of opsins. Some TGD paralogs, including those of parapinopsin (PP), Vertebrate Ancient (VA) opsin, and long wavelength-sensitive (LWS) opsin, show different absorption spectra and/or expression patterns. However, our knowledge of detailed evolutionary processes and mechanisms by which TGD contributed to opsin diversification is still limited. Here, we report that TGD paralogs of parietopsin (PT1 and PT2) are retained in extant species. While most species, including the zebrafish, have only PT1, several characins, including the red piranha (Pygocentrus nattereri), retain both PT1 and PT2. The Mexican tetra (Astyanax mexicanus) and catfishes have only PT2. To assess the degree of functional diversification between PT1 and PT2, we characterized spectral properties and expression patterns. Maximum absorption spectra differ slightly among PTs. Those of red piranha PT1, PT2, Mexican tetra PT2, and Japanese catfish (Silurus asotus) PT2 were located at 517 nm, 528 nm, 517 nm, and ~535 nm, respectively. Fluorescence in situ hybridization showed that (1) piranha PT1 and PT2 are expressed in the same pineal cells, and (2) they are also co-expressed with PP1. Based on these findings, functional diversification of opsins after gene duplication is discussed.

Universal Single-Copy Orthologs (USCOs) are taxonomically-defined gene sets originally proposed to evaluate completeness in genome assemblies. Due to their single- or low-copy nature, they have also been proposed as ideal markers for phylogenomics. Here, we test this assumption using two vertebrate datasets, mammals and ray-finned fish, and their corresponding USCO reference marker sets, actinopterygii_odb10 and mammalia_odb10. Our results show that although most USCO markers are single-copy, 85% and 91% respectively, empirical analyses using a large number of taxa are able to identify a non-negligible proportion of paralogs in the actinopterygii_odb10 and mammalia_odb10 BUSCO gene sets. Gene phylogenies show that several of these paralogs derive from ancient whole-genome duplication events, particularly the teleost-specific whole-genome duplication, and also more recent duplications such as those in salmonids and cyprinids. Paralogs do not reflect the organismal phylogeny but rather the history of gene duplications, so retaining paralogs in phylogenomic datasets can lead to inaccurate phylogenomic hypotheses. To test the effect of ignoring paralogy in BUSCO gene sets, we analyze the mammal and ray-finned fish datasets using concatenated maximum likelihood and summary coalescent methods by randomly choosing one paralog per locus if these are present. The resulting phylogenomic trees are overall consistent with the reference phylogenies, although small discrepancies were found. Phylogeny-guided paralogy removal also changes the sequences retained for many loci, but whether these changes alter the inferred species trees remains to be tested.

Parasitoid wasps develop in or on other arthropods and have evolved a wide range of strategies for exploiting their hosts. Ascogaster reticulata is a solitary egg to larval endoparasitoid wasp belonging to the subfamily Cheloninae of Braconidae. We constructed a genome assembly using long-read sequence data. The assembled genome was 280.5 Mbp, consisting of 209 contigs. The N50 value and GC content ratio were 13.9 Mbp and 30.43%, respectively. The percentage of transposable elements and repetitive regions was 52.43%. The predicted number of genes is 27,261. Merqury revealed that the QV value and completeness of the genome were 61.79 and 99.3917%, respectively. BUSCO analysis detected approximately 99% of the core genes. These results suggest that the genome data were of sufficient quality for use as a reference genome. Therefore, these data can be utilized as a genomic resource for comparative studies within Braconidae and for future investigations into the biology and evolution of Cheloninae parasitoid wasps.

This study introduces a novel model that integrates the capacity of Potts models to capture structural dependencies with the computational efficiency afforded by Felsenstein's algorithm in phylogenetics. While the present work focuses on protein models, the PhyloPotts framework is equally applicable to alignments of discrete states, including RNA, DNA, and 3DI characters. Our findings indicate that, although phylogenetic estimates do not consistently exhibit reduced bias relative to standard models, the approach substantially enhances the accuracy of ancestral sequence reconstruction. The method is freely available via an open source implementation in the PhyloPotts package for BEAST 2 from https://github.com/rbouckaert/potts.

GC-biased gene conversion (gBGC) is a widespread recombination-associated process that favours the fixation of GC alleles independently of their fitness. In coding regions, it can oppose purifying selection and favour the accumulation of deleterious AT to GC mutations. In mammals where recombination hotspots are predominantly determined by PRDM9, the deleterious effect of gBGC is expected to be episodic due to the rapid turnover of hotspots. The deleterious impact of gBGC episodes on proteins has been well documented, particularly in rodents and primates. However, the evolutionary response of the affected proteins after the episodes end has been scarcely investigated. Theory predicts that a decline in gBGC intensity should trigger a compensatory response that enables previous damage to be repaired. Recently, we reported evidence for pervasive compensatory evolution following gBGC episodes in Hydromyini rodents. Here, we aimed to investigate whether this post- gBGC compensation is widespread among mammals in which gBGC is episodic. Using a phylogenetic approach, we analysed substitution patterns in coding exons from five mammalian taxa spanning a range of effective population size (Ne): Murini, Rattini, Cetacea, Cercopithecidae and Hominoidae. We detected thousands of gBGC episodes, characterised by an increased rate of non-synonymous AT to GC substitutions. We found evidence for compensatory evolution in non-primate taxa, and a lack of signal in Primates. Although the influence of Ne remains unclear, our study highlights both the evolutionary significance of post-gBGC compensation and the methodological challenges in detecting it.

Whether the homologous parallel triplex (HP-triplex, formed by dsDNA and a homologous ssDNA) acts as an intermediate prior to DNA recombination remains controversy, owing to the intrinsic instability and lake of experimental data. Here, we report a HP-triplex that is stable at least at 35 oC under protein-free conditions. Two complementary circular ssDNAs were first hybridized at a 1:1 ratio to form an LR-chimera containing both Z-DNA and B-DNA. Then the circular homologous ssDNA (T-rich) was found to easily hybridize with this LR-chimera to form the HP-triplex (strand ratio of 2:1), involving both left- and right-handed triplex parts under the topological constraint. The third strand locates in the major groove and forms hydrogen bonds with both bases in each base pair of dsDNA without obvious strand exchange. Our data indicate that the third-strand ssDNA may first bind to the B-DNA region, followed by association with the Z-DNA region to assemble the full triplex structure. Structural simulations and analyses reveal that topological constraint imposed by the circular architecture raises the energy barrier for strand displacement, thereby stabilizing this triplex. In summary, the dynamic HP triplex can be stabilized via a Z B junction structure that mimics the in vivo one formed by negative supercoiling. These findings demonstrate that HP-triplex may act as an intermediate for DNA recombination, anti-gene strategy, and other interactions between dsDNA and homologous ssDNA or RNA.

Biological populations have been evolving in multi-dimensional trait spaces. Ecological interactions can induce the evolutionary diversification of a biological population into distinct multiple populations, which is called evolutionary branching in adaptive dynamics theory. If possible mutations are restricted to certain directions due to developmental, physiological, or physical constraints, the resulting evolutionary dynamics may be restricted to subspaces with fewer dimensions than the original trait spaces. In general, such constraints are also evolving slowly, due to rare and slight mutational deviations from the constraint subspaces. In this sense, constraints are not strict. However, the method for describing non-strict constraints that is globally consistent has yet to be developed. This study proposes a general and flexible method for constructing a non-strict constraint subspace from the corresponding strict constraint subspace. This method enables efficient analysis of adaptive evolution including evolutionary branching under non-strict constraints in the framework of adaptive dynamics theory. Applying this method to a resource competition model shows that the shape of the carrying capacity distribution can be treated as a non-strict constraint curve, and that the likelihood of evolutionary branching is affected not only by the shape of the constraint curve but also by its slow evolution.

The evolution of genome size and eukaryotic genome organisation has puzzled evolutionary biologists for decades. Are differences in DNA content and protein-coding capacity of eukaryotic genomes a result of natural selection, or neutral processes? If they were shaped by natural selection, which features might be connected with these differences: cell volume, ecological strategy, availability of nutrients, etc.? Kinetoplastids are protists with compact nuclear genomes that almost entirely lack introns, but whose mitochondrial DNA content (kDNA) can rival that of their nuclei. We hypothesised that the reduction in non-protein coding nuclear DNA content during the evolution of kinetoplastids is connected with an increase in kDNA content: species with larger fractions of total cellular DNA that is allocated to kDNA will have concomitantly smaller fractions of non-coding DNA in the nuclear genome. We developed a framework for estimation of DNA content in both the mitochondria and nuclei from quantitative fluorescence microscopy data, and inferred the fraction of non-coding DNA in nuclear genomes based on the combination of transcriptomic and fluorescence imaging data. Across 10 genetically distant isolates from diverse free-living kinetoplastid lineages, there is indeed a trend towards a lower proportion of non-protein-coding nuclear DNA when kDNA proportion is higher relative to the total DNA in the cell. This observation is consistent with the hypothesis that in free-living kinetoplastids, mitochondrial and nuclear genomes might co-evolve in terms of their DNA content. This finding raises the possibility that natural selection acts on genome size evolution in kinetoplastids, and contributes evidence and methodology to address this question more broadly in eukaryotes.

Emergence and evolution of functional RNA and protein structures are the central problems for understanding the origin of life. Although it is well known that catalytically active RNA elements, ribozymes, can catalyze many reactions, including peptide bond formation, the specifics of the transition from the hypothetical, primordial RNA world to protein-based life centered at the translation system remain enigmatic. We developed AMES, Atomistic Molecular Evolution Simulator, and employed it to perform computer simulations of the evolution of short RNA molecules and RNA-peptide complexes. Comparison of the evolutionary trajectories and the structures of RNA molecules and RNA-peptide complexes emerging in these simulations shows that short random peptides accelerate RNA evolution, stabilize RNA folds, and boost the structural diversity of evolving RNA molecules, potentially enabling a broader range of activities. We hypothesize that the RNA world was, actually, an RNA-peptide world, in which, from the earliest stages of evolution, evolving RNA molecules interacted with short random peptides synthesized in a non-templated manner. These interactions could drive the evolution of diverse RNA structures and activities, and eventually, of the translation machinery.

Hox proteins are evolutionary conserved transcription factors guiding morphogenesis in all bilaterian animals. Conservation mainly concerns the homeodomain (HD), a helix-turn-helix DNA binding motif and flanking sequences close by. Besides, sequences are poorly conserved and disordered. Here we investigate how poorly conserved Intrinsically Disordered Regions (IDRs) influence Hox protein activity, focusing on a closely related Drosophila Hox pair and using a wide range of in vivo function. Studying Drosophila Ultrabithorax (Ubx) and AbdominalA (AbdA) proteins, we find that conserved protein motifs close to the HD are differently used in Ubx and AbdA. Coarse-grained molecular dynamics simulations identify a large set of intra-molecular interactions, notably involving poorly conserved disordered regions upstream and downstream the HD. Swapping these regions in Ubx and AbdA alters protein functions and switches motif dependency in some biological situations. This suggests that IDRs, predicted to mediate large and dynamic sets of intramolecular interactions, influence protein motif usage in a context dependent manner.

Forward-in-time population genetic simulations are widely used in evolutionary analyses, but simulating large populations and long genomic regions remains computationally demanding. To reduce this cost, parameter rescaling is widely employed, in which the original evolutionary process is approximated by one with a smaller population size and fewer generations. Recently, several studies using the SLiM simulator have raised concerns about the accuracy of this rescaling approach. In this study, we show that many of the biases reported in these studies can be mitigated by using a different simulation algorithm. These results reveal that the accuracy of parameter rescaling depends on how well the simulation algorithm preserves diffusion-limit properties under rescaling.

Aquaporin-10 (AQP10) is an aquaglyceroporin that transports small, uncharged molecules, such as water and glycerol. In some species, such as mice and cattle, the AQP10 gene (Aqp10) is pseudogenized. In rodents, pseudogenization of Aqp10 is limited to the Myomorpha suborder, while Aqp10 remains intact in other rodent lineages. The order Cetartiodactyla encompasses diverse species, including Suina, Tylopoda, Ruminantia, and Whippomorpha. The distribution of intact Aqp10 across the entire Cetartiodactyla and analyses estimating the timing of Aqp10 pseudogenization have not been sufficiently investigated. In this study, we identified species within Cetartiodactyla that retain intact Aqp10 and analyzed the detailed patterns of pseudogenization to elucidate the specifics of pseudogenization within this order. We collected Aqp10 loci by BLAST and synteny analyses of publicly available genome databases of 51 cetartiodactyl species from 17 families. Dot plot and alignment analyses revealed that seudogenization of Aqp10 occurred in 42 Cetartiodactyla species, including the Camelidae, warthog (Suidae), Antilocapridae, okapi (Giraffidae), Cervidae, Bovidae, Hippopotamidae, and Cetacea. In contrast, nine species retained an intact Aqp10: pig and babirusa (Suidae), peccary (Tayassuidae), mouse deeres (Tragulidae), giraffes (Giraffidae), and musk deers (Moschidae). The timing of pseudogenization was estimated by comparing patterns of exon deletions, frameshifts, and nonsense mutations in Aqp10 pseudogenes. RT-PCR analysis of cattle and sheep tissues revealed that these species express Aqp3 and Aqp7 in their intestines, but not Aqp10. Collectively, these findings indicate that Aqp10 pseudogenization occurred independently in multiple cetartiodactyl lineages and provide new insights into the evolutionary history of Aqp10 loss within this order.

The allorecognition system of the hydroid Hydractinia symbiolongicarpus plays an essential role in mediating self vs. non-self recognition among cnidarian colonies, analogous to the role of the major histocompatibility complex (MHC) in bilaterian immune systems. Both complexes have contiguous stretches of similar, highly duplicated genes, suggesting a homologous relationship. Protein structural prediction methods such as AlphaFold2 also indicate structural similarity between the proteins encoded within the Hydractinia ARC and the MHC. We have determined the structure of the ARC in the related hydrozoan Podocoryna americana to determine the degree of conservation of the cnidarian allorecognition complex. Protein structural prediction methods have revealed a number of candidate allorecognition proteins that form a long contiguous complex. Phylogenies suggest homology between cnidarian allorecognition genes, as well as syntenic conservation of this gene complex across hydrozoans. Ongoing comparative genomic analyses across bilaterian species aim to clarify the relationship between the ARC and the MHC, thereby informing our understanding of the evolutionary origins of bilaterian adaptive immunity.

MicroRNAs have been likened to the "dark matter" of eukaryotic genomes, reflecting their pervasive regulatory influence. MicroRNAs were first identified through genetic studies of developmental timing in the nematode Caenorhabditis elegans. Let-7 was the first microRNA recognized to be broadly conserved. The principal target of Let-7 in the developmental timing pathway is the TRIM-NHL RNA-binding protein LIN-41. During the L4 larval stage, Let-7 represses lin-41 translation by binding to two Let-7 complementary sites in the lin-41 3'UTR. Despite the importance of microRNA-based translational regulation, the underlying molecular mechanisms are incompletely understood. Through genetic analysis, we discovered an unrecognized feature of the mechanism by which Let-7 controls lin-41 translation. This mechanism requires a 5'-regulatory exon containing a seven-amino acid upstream open reading frame (uORF) and conserved sequence elements. Genome editing indicates that the specific uORF amino acid sequence itself is not important. Our data suggest that uORF translation and 5'UTR structure limit initiation at the downstream lin-41 start codon, enabling tight control by Let-7. Without this mechanism, the Let-7 microRNA is unable to properly regulate lin-41 to enable proper development.